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Abstract 



Symplectic ensemble of disordered non-Hermitian Hamiltonians is studied. 
Starting from a model with an imaginary magnetic field, we derive a proper 
supermatrix fi-model. The zero-dimensional version of this model corresponds 
to a symplectic ensemble of weakly non-Hermitian matrices. We derive ana- 
lytically an explicit expression for the density of complex eigenvalues. This 
function proves to differ qualitatively from those known for the unitary and 
orthogonal ensembles. In contrast to these cases, a depletion of eigenvalues 
near the real axis occurs. The result about the depletion is in agreement with 
a previous numerical study performed for QCD models. 
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I. INTRODUCTION 



Properties of non- Hermit ian random operators or matrices have attracted recently a 
considerable attention. Non-Hermitian random Hamiltonians can appear as a result of 
mapping of a model for flux lines in a (rf + l)-dimensional superconductor with line defects 
in a tilted magnetic field on a (i-dimensional model for bosons in a random potentiali. Non- 
Hermitian operators enter Fokker-Planck equations that describe diffusion and advection 
of classical particles in a spatially random but time-independent velocity fieldH and also 
determine equations used for study of problems of turbulenceOcl. 

Ensembles of random complex non-Hermitian and real asymmetric matrices find their 
application for a description of dissipative quantum maps0 in neural network dynamics^. 
Recently it was suggested that they could be relevant for QCD where they correspond 
to a random Dirac operator with a non-zero chemical potentialS Starting from the first 
workslll0, properties of the ensembles of the non-Hermitian matrices were intensively studied 
in a considerable number of publicationJli^. 

Unusual properties of the ensembles of the non-Hermitian operators or matrices are re- 
lated to the fact that eigenvalues of the operators and matrices can be complex. Completely 
different methods have been used for study of distributions of the eigenvalues on the com- 
plex plane. For example, the authors of Refs.Uli applied Green functions methods while in 
Refs.i2l, a method of orthogonal polynomials was used. In Ref.El, a new regime of a "weak 
non- Hermit icity" was found and the authors have calculated a joint probability of complex 
eigenvalues for complex weakly non-Hermitian matrices. For calculations they used the su- 
persymmetry technique^ and derived a zero-dimensional non-linear cx-model. An important 
information about the distribution function of complex eigenvalues oi N x N matrices for 
the orthogonal, unitary and symplectic chiral random matrix ensembles has been obtained 
recently numericallyil. 

Although the model with the non-Hermitian Hamiltonian of Ref.S differ from those with 
random matrices, they turn out to be closely related to each other. In Ref.ii, the model 
with the non-Hermitian Hamiltonian H was studied using the supersymmetry method. This 
Hamiltonian can be written in the form 

,1.1, 

where p = —iV, m is the mass of particles, and U (r) is a random potential. The vector h 
is proportional to the component of the magnetic field perpendicular to the direction of the 
line defects in the initial problem of the vortices in superconductors. The Hamiltonian ( |1 . 1|) 
describes a particle moving in an imaginary vector-potential ih and a real random potential 
U (r). The distribution of complex eigenvalues on the complex plane can be extracted from 
the distribution function P(e, y) defined as follows 



Pie,y) = Tj{y:6{e-e',)5{y-el)) , (1.2) 




where e'/^ and e'l are the real and imaginary parts of the eigenenergies, respectively and V 
is the volume of the system. The angle brackets stand for an averaging over the random 
potential and the sum should be taken over all states. 
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The problem of calculation of the function P(e, y) was mapped in Ref.t3 onto a new 
supermatrix non-linear cr-model. This model differs from the conventional ones written 
previously0 by the presence of new "effective fields" . The symmetry of the matrix Q entering 
the new cr-model is the same as that of obtained in Ref.0 for the orthogonal ensemble. This 
is not accidental because the Hamiltonian, Eq. ( |1 . 1| ) , is real and, hence, time reversal 
invariant. 

To violate the time reversal invariance one can add to the Hamiltonian ( |1 . ID a real 
magnetic field or/and magnetic impurities. This leads to additional terms in the a-model 
lowering the symmetry of the model. As a result, one gets@ the a-model with the super- 
matrices Q corresponding to the unitary ensemble. Although the real magnetic interactions 
in the Hamiltonian with the imaginary vector potential do not correspond to any physical 
interactions in the initial problem of the vortices, consideration of the cr-model for the uni- 
tary ensemble was interesting from the formal point of view because it allowed to establish 
important relations with the random matrix models. 

The cr-models corresponding to the Hamiltonian of Eq. (|1 . 1| ) and its extensions can be 
written in an arbitrary dimension. Remarkably, the zero- dimensional version of the a-model 
for the unitary ensemble is exactly the same as the zero-dimensional cr-model derived in 
Ref.Ei for the weakly non-Hermitian matrices. Complex random non-Hermitian matrices 
appeared in studies of dissipative quantum map^^lll, which justifies an interest to studying 
the unitary ensemble. By the term "weakly non-Hermitian" the authors of Ref.El called 
matrices X that could be represented in the form 

X = A + iaN-^^^B , (1.3) 

where A and B are Hermitian N x N matrices and a is a parameter characterizing the 
non-Hermiticity. 

The cr-model obtained from the ensemble of the matrices, Eq. (|1.3D, and from the 



Hamiltonian with the imaginary vector-potential allows us to relate the parameters h and 
a to each other. A similar correspondence exists for the orthogonal ensemble. 

Study of the distributions of the complex eigenvalues revealed a striking difference be- 
tween the orthogonal and unitary ensembles. The function P{e,y), Eq. (|1.2| ), is a smooth 
positive function of y for the unitary (provided the disorder is not very strong, this function 
does not depend on e). It reaches its maximum at y = and monotonously decays with 
increasing y. The corresponding function P(e, y) for the orthogonal ensemble is a sum of a 
smooth function and a (5-function of y. This means that a finite fraction of the eigenvalues 
remain real at any degree of the non-Hermiticity. 

In all the works done in statistical physics only the orthogonal and unitary were con- 
sidered. The symplectic ensembles have not been even mentioned, apparently due to the 
absence of any applications. However, for the random matrix ensembles applied for clarify- 
ing properties of QCD model the symplectic ensemble is of the same importance 
as the orthogonal and unitarvones. Moreover, numerical results for distributions of complex 
eigenvalues presented in Ref.Ba demonstrate a pronounced difference between the ensembles. 
The distribution of the complex eigenvalues on the complex plane is homogeneous in the case 
of the unitary ensemble while it shows an accumulation of the eigenvalues along the real axis. 
This corresponds to the presence of the 5-function in the function P{e,y), Eq. (|I.2[), found 



3 



in Ref.0. Although the authors of Ref.il considered chiral matrices, the dependence of the 
number of eigenvalues at the real axis on a parameter characterizing the non-Hermiticity 
was found to be exactly the same as in Ref.S. This shows that the phenomenon of the 
accumulation is quite general. 

A completely different behavior was found in Ref.ii for the symplectic ensemble. The 
distribution function of the complex eigenvalues is in this case smooth but the probability 
of real eigenvalues turns to zero, which corresponds to a depletion of the eigenvalues along 
the real axis. This is a new effect that clearly motivates an analytical investigation of 
non-Hermitian symplectic matrices. 

In the present publication the distribution function P{e,y), Eq. ( |1.2| ), is calculated for 
the ensemble of symplectic non-Hermitian matrices. This is done by writing a proper zero- 
dimensional (T-model. We are able to obtain an explicit expression for the function P(e, y) 
and demonstrate the depletion of the eigenvalues along the real axis. 

The paper is organized as follows: In Sec. II, we introduce the notations and remind 
to the reader the scheme of the derivation of the cr-model. In Sec. Ill, we present the 
parametrization of the supermatrices Q for the symplectic ensemble. In Sec. IV, the joint 
probability density of complex eigenvalues is calculated. Sec. V summarizes the results, and 
in the Appendix the Jacobian of the parametrization is derived. 



II. NON-LINEAR cr-MODEL 

The derivation of the a model for the non-Hermitian orthogonal and unitary ensembles 
has been comprehensively presented in Ref.@. Addressing to this paper for all details, 
we repeat some intermediate steps, concentrating on minor changes that have to be done 
in the symplectic case. The final goal is to derive the joint probability density of complex 
eigenenergies P(e, y), Eq. ( |1.2| ). Of course, one can derive the zero- dimensional cr-model from 
the ensemble of symplectic random matrices but we prefer to start from the Hamiltonian, 
Eq. (|1.1|) , adding to it spin-orbit impurities. 

Due to the non-Hermiticity the Hamiltonian, the notion of advanced and retarded Green 
functions, and Gf usually used in perturbation theory and in deriving the non-linear 
cr-models becomes meaningless since they loose their analytic properties. The difficulty can 
be overcome by introducing an Hermitian double size operator M of the form 

M-( Z7'.^.r^Zr'l], (2.1) 



where 



-i{H"-y) -(H'-t) 



H'=^"^"*\ H'. ^ , (2.2) 



In equations ( |2.1| ) and ( |2.2| ), H is the Hamiltonian, Eq. ( |1.1| ), and means its Her- 
mitian conjugated. Instead of manipulating the non-Hermitian operator, one can use the 
Hermitian operator Using the "effective Hamiltonian" M, Eq. ( |2.1| ), one can represent 
the complex eigenvalues distribution function P{e,y), Eq. ( [L.2|) , in a form of a functional 
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integral over supervectors ip (r) with the weight exp(— £) with the Lagrangian C taking the 
form 

C = -tJ i^{r)[no + U{r) + Ko (r)]^(r) dr . (2.3) 

Here, '?/'(r) and ip(r) are the standard supervector and its charge-conjugated counterpart, 
respectively, composed from anticommuting and commuting fields0. The matrix operator 
Ho consists of two terms 

no = {K-e + z^A)I + iK^{H'i + yra) , (2.4) 

where H'q and Hq have the form 

In equation ( |2.4|) , 7 is a small positive number that should be put to zero at the end of 
calculations. The term Vso in Eq. (p.3| ) stands for spin-orbit impurities. It can be derived 
from the initial Hamiltonian after a formal inclusion of the interaction Ugo (i") with the 
spin-orbit impurities. The simplest form of this interaction can be written as follows 

Uso{r) = CT [Vn3o(r) x p] , (2.6) 

where the vector <t is formed from the Dirac matrices a^-, c^, and a^- The matrices J, A, Ai 
and Ts entering Eq. ( p.4|) have the form 

'-[\\\ ^AVX ^-[\-X 




" ■ (2.8) 



Due to the necessity of considering the spin variables the supervectors ^/^ (r) have now 16 
components. The unit blocks 1 in the matrices in Eq. (|2.7| ) have the size 8x8 and unities 
1 entering the matrix ra are 2x2 matrices. 

The distribution of the electric fields Vuso(r) is assumed to be Gaussian: 

(Vwso(r)) = , (9,Wso(r) 9,w,o(r)) = '^^^'f ~ ""^^ , (2.9) 

where the density of states at h = at the Fermi surface v = mpo/27r^ and is the 
spin-orbit scattering time. 

Further transformations are performed according to standard rules of the supersymmetry 
techniqueiil. Averaging over the disorder results in an interaction term ip'^ in the Lagrangian 
C This term is decoupled by integration over a supermatrix Q. Then, one integrates the 
supervector i{j out, arriving thus at an integral over Q with the weight exp{—F[Q]), where 
F[Q] is a free energy functional. 
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The spin-orbit interactions lead to additional "effective fields" of a certain symmetry in 
the free energy functional These fields lower the symmetry of the functional. As a 

result, a part of fluctuational modes have a gap and their contribution at low energies can 
be neglected. This is equivalent to putting certain elements of the supermatrix Q to zero. 
Carrying out this procedure, one comes to a matrix Q with spin blocks proportional to unit 
matrices. This is equivalent to consideration the model with 8x8 supermatrices Q having 
a new symmetry. These are the same supermatrices as those used in Ref.0 for description 
of the symplectic case. 

Therefore, we should perform calculations similar to those of Ref.@ but integrating over 
the supermatrices Q with the symmetry corresponding to the symplectic ensemble. Af- 
ter standard transformations, we reduce the distribution function P(e, y) to the following 
integral 

^(e, y) = -^J ^[Q] i-F[Q]) dQ ' (2-10) 



"^[Q] — ('542 + ^542) ('524 + '^24) ~ ('^42 + '542) ('^42 + Q 

with the zero- dimensional version of the free-energy functional 

F[Q] = STr (^[Q, Ai]' " f AiTsq) • (2.11) 

In equation ( |2.11|) , the symbol [..] stands for commutator, STr for supertrace and we have 
introduced the following parameters: 

o^^^, .^f^, (2.12, 

where Dq is the classical diffusion coefficient and A = {2vV)^^ is the mean level spacing 
(the factor 2 in this expression is due to hfting of the spin degeneracy by the spin-orbit 
impurities) . 



III. PARAMETRIZATION FOR THE SUPERMATRICES Q 

To integrate over all symplectic matrices, proper variables parametrizing the Q super- 
matrix should be introduced. The parameters have to be chosen so that to cover all the set 
of the symplectic matrices: Any symplectic matrix has to be reached only once. 

Although the parametrization for the unitary and the orthogonal ensembles cannot be 
used for the symplectic ensemble only minor changes have to be done to adjust the non- 
Hermitian parametrizationB to the case under consideration. As in Ref.il, we represent the 
Q matrix in the form of the product 

Q = TYQoYT . (3.1) 

To fulfill the constrain = 1, the following equalities must be hold 
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Ql = l , TT = 1 , YY = 1 . (3.2) 

As in Ref.ii, the supermatrices T and Y are chosen to commute with Ai 

[T,Ai] = [r,Ai] = 0. (3.3) 

The next simphfication facihtates the parametrization and enables one to calculate the 
Jacobians quickly. Namely, we decompose the supermatrix Y into the product of the matrix 
Yq containing commuting variables and the matrices R and S, consisting of the Grassmann 
ones 

Y = YqRS . (3.4) 

The 2x2 blocks in the matrices R, S, and Qq are chosen to be diagonal; the necessary 
symmetry of the 2x2 blocks a, b and a is achieved by a proper choice of 2 x 2 blocks of the 
matrix Yq. Thus, the matrices Qq, R, and 5* can be written in a form similar to the one inll 

y—Tssimp —cos if J \0 ix J 

where ip and x proportional to the unity 2x2 matrix, 

R^f^O], fi-fl-^'-P ). (3.6) 

[or \ -2p l + 2ppr ^ ' 



and 



2^a l-2aO' "=Uo " ^^'^^ 



The matrices p and a in Eqs. ( p.6| , p.7| ) have the form 



The next step is to represent the supermatrix Yq in Eq. ( |3.4| ) as the product 

Yq = Y^Y^Y^ , (3.9) 

where Y^ is the diagonal matrix 

/exp(^/3/2) \ A_f/5r3 \ 

exp(z/3/2)j' ^-l, Arsj- ^^'^^^ 

In order to recover the symplectic symmetry we have to choose the matrices Yi and Y2 
as follows 













, w = 







w 1 




w 1 



n .. I ^ = n > (3-11) 



w 



cosh(/i/2) — i sinh(/i/2) \ 
zsinh(/x/2) cosh(/x/2) J 



and 



^ / cos(^2/2) -^sm{9,/2)\ 
' \-ism{92/2) cos(^2/2) / ' ^ ' ^ 



^'-[ oj' ^'~[io 

Finally, the supermatrix T can be taken as 

u 



, . W cos(^/2) sin(^/2) \ / \ 
' l^-zsin(^/2) cos(^/2) / ^ / ' 



where 

e 



e 

i^i 



fl-2r]r] 2r] \ fl-2KK 2k \ 

-2^ 1-2^W' -2^ 1-2^.J- ^3.14) 

The 2x2 matrices 6 and in Eq. ( p.l3|) are proportional to the unit matrix and the 
matrices t] and k, are 

The explicit form of the supermatrix Q within the parametrization suggested, Eqs. (|3.5| - 
3.14| ), is very similar to that for the orthogonal ensemble@ and differs from the latter by 
minor changes in the matrices w, 62, cr, p, rj, and n. 

To ensure the unambiguity of the parametrization, we should specify the variation range 
of the variables. This is done by comparing the compact and the noncompact sector with 
those in the standard parametrizationEll. As the result, the variables vary in the following 
intervals: 

- 7r/2 < 99 < 7r/2 , 0<x<oo, -n < 9 < n , -00 < 6*1 < 00 , (3.15) 
< jj, < TT , 0<6'2<oo, < (3 < IT , < /3i <2tt . 

The only thing that remains to be done to perform explicit calculations for physical 
quantities is calculation of the Jacobian of the transformation to the variables described by 
Eqs. 

Its derivation presented in the Appendix leads to the following final result for the ele- 
mentary volume 

[dQ] = J^JgJ^JcdRB dRp , (3.16) 

where 
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_ 1 cosv^coshx ,1 1 



Stt (sinhx + ^sinyj)^ ' 327r sinh ^{9i 



lb 



J 1 sin 62 sinh /i ^ 4 sinh^ x 

2^7r2 (cos 6*2 — cosh /i)^ ' (sinh x~ sin (^9)2 

and 

(IRb = dO dOidif dxd^ dd2d(3 d(3i , = drj drj* dn dK* da da* dp dp* . (3.18) 

Equations ( p.l| - p.l7| ) are sufficient for evaluation of any integral over the supermatrices 
Q and, with the help of Eqs. ( p. 10 , 2.11 ), provides a straightforward way of calculating the 



distribution function of complex eigenvalues P{e,y), Eq. (JT^. 

IV. DENSITY OF COMPLEX EIGENVALUES 

Before staring the calculations let us introduce more compact notations. As it will 
be seen in what follows, only the following combinations of the variables describing the 
parametrization of the supermatrix Q enter all functions of interest 

t = simp , z = sinh x , 1^ = cosh p , X = cos 62 . (4.1) 

Since the matrices T and Y commute with Ai, the ffist term in the free energy, Eq. 
( p. lip does not depend on them. The second term in Eq. ( p.ll| ) does not depend on T. As 



a result, the free energy F[Q] takes a rather simple form 

F[Q] = a^{f + z^) + x\{\t - iuz) + A{aa* + pp*){uj - \){t - iz)] . (4.2) 



The fact that F[Q\ does not depend on T simplifies the integration over Q in Eq. ( p. 101) 



Using the parametrization, Eqs. ( |3.1|3.13| ), we can also represent the supermatrix Q as 

Q = uQu , (4.3) 



with u from Eq. ( |3.14| ) and some supermatrix Q. Substituting Eq. (|4.3| ) into Eq. ( |2.10|) for 



the density of complex eigenvalues P(e, y) and integrating over 77 and r/*, one represents this 
function in the form 



P{^. y) = ^J [Str{T,A^Q)]' exp(-F[Q]) dQ (4.4) 
exp(-F[Q]) dQ . 



Anu d^ 
^ A dx^ 

For the symplectic ensemble, one has in Eq. (^4.4| ) an uncertainty of the type Ox 00, since 
the integrand does not contain the variables k and n* and, on the other hand, the Jacobians 
Jq and are singular as 6*, 6*1, 6*2, and /i — 0. To resolve this singularity, we can use the 
regularization procedure developed for the orthogonal ensemble". All the manipulations are 
identical to those of Ref.0, because the free energy, Eq. (|4.2| ) has the form similar to that of 
the orthogonal ensemble. Moreover, the singularities of the Jacobians, Eqs. ( p.l7| ), are the 
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same as the ones for the orthogonal ensemble. We do not specify the procedure once more 
and present here only the final result of the regularization with proper changes of notations. 
The function P(e, y) can be written in the form of a sum of two terms 

P(e,y) = P«(e,y) + p(2)(6,y), (4.5) 

where 



and 

^^'ne,l/) = ^^/exp[-a2(t2 + z2)-x(tu;-az)] (4.7) 

it - izY z^x^ , 7 , 
X — TT, TTT^ — at az duj aX . 

(t2 + ^2)2 

The integration in Eqs. ( [4.6|) and ( [4.7|) is performed in the intervals — 1 < t < 1, — oo < z < 
oo, 1 < < oo, and — 1 < A < 1. 

To perform the integration in Eq. ( |4.7] ) over A, one should introduce an infinitesimal 
positive 6, defining z_ according to Z- = z + i6sgn{x), so that the integral becomes conver- 
gent. After this, the integration over A and uj in Eq. ([4.7|) is easily carried out. Adding the 
result of the integration to Eq. (^I6|) we obtain 

ly r+1 r+oo 

p(^e,y) = — — J dt J dzexp{-a^{t'^ + z'i))[{t + iz^fexp{ixz^-tx) (4.8) 
— {iz_ — tY exp(zxz_ + tx)] 



2 \2 



it{t^ + z: 

Comparing Eq. ( [4.8| ) with its analog for the orthogonal ensemble, we notice an important 
difference between them: The variable Z- is present in the numerator in Eq. ( ^4.8| ), whereas 
it stands in the denominator of the equation for the orthogonal ensemble. In the latter 
case, the distribution function P{e,y) contains an additional contribution of a 5-function 
after differentiation of z_ over x. For the symplectic ensemble, the differentiation of z^ in 
the numerator leads to no singularity on the real axis and one take the limit 5^0 before 
calculating the integral in Eq. (|4.8| ). Thus, only the exponents should be differentiated over 
X. These differentiations simplifies considerably the integrand and the integration over z 
can be easily carried out. After that one obtains 

, u X ^ f \ r / 9 9x sinh(tx) 
^(^■^) = AS?^^^n-i?U ^'-M-aV)^. (4.9) 



Equation (|4.9| ) solves completely the problem involved and is the main result of the 
present paper. 

The following properties of the density function P(e, y) are easily checked: It is symmetric 
with respect to y and is properly normalized 
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J dyP{e,y) = l. (4.10) 

In the limit a ^ l(the limit of a strong non-Hermiticity), one obtains the universal 
asymptotics valid for all three ensembles 

The form of the density of complex eigenstates, Eq. ( [4.11|) , corresponds^ to the "elliptic 
law" of Refs.ffl. 

In the opposite limit a <^ 1, the function P{e,y) can be written as 



2 / 2 

U X ^ IX 



^(v!/) = ^J^v/^exp -j^ . (4,12) 



The behavior of the function P(e, y), Eq. ([4.9|), at small y (related to x by Eq. (|2.12|)) is 



drastically different from the behavior of the corresponding functions for the orthogonal and 
unitary ensemblesii. This function is small at small y, being proportional to y'^ and turns 
to zero in the limit y ^ 0. This means that the probability that eigenvalues remain real 
at finite degree of the non-Hermiticity is zero. In other words, the distribution function of 
complex eigenvalues exhibits a depletion along the real axis. The depletion region broadens 
with increasing the non-Hermiticity. The function P(e, y) is represented in Fig. 1 for several 
values of a ~ 1. 



V. CONCLUSIONS 

In the present paper, we studied analytically disordered non-Hermitian models with 
the symplectic symmetry. This is the last of three universality classes that has not been 
considered yet. Using the supersymmetry technique we derived a proper non-linear a- model 
starting from the a model of disorder with a direction. The zero-dimensional version of the 
non-linear cx-model corresponds the ensemble of random non-Hermitian symplectic matrices. 
Within the zero-dimensional cx-model, we calculated the joint probability density function 
of complex eigenvalues. We introduced a convenient parametrization and calculated the 
Jacobian corresponding to this parametrization. 

All this allowed us to derive an explicit expression for the density of complex eigenvalues. 
Asymptotic behavior of this function demonstrates clearly that the basic properties of the 
system depend strongly on the ensemble. Introducing the non-Hermiticity in the Hamilto- 
nian affects very differently the spectrum of three ensembles. Only when the non-Hermiticity 
is very large, the difference is no longer important. 

It is known from previous works that the eigenvalues of a system belonging to the uni- 
tary ensemble are smoothly distributed around the real axis. The density function for the 
orthogonal ensemble contains a 5-function contribution on the real axis describing an accu- 
mulation of the eigenvalues. In contrast to the previous cases, we obtained for the symplectic 
ensemble a depletion of eigenvalues along the real axis, which is in a good agreement with 
the results of a numerical studyc3. These features correspond to a tendency of a system 
from the orthogonal ensemble to preserve localized behavior. However, after introducing 
spin-orbit impurities, the system acquires delocalized features. 
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VI. APPENDIX 



The Jacobian of the parametrization specified by Eqs. ( |3.1[ - ^.14] ) can be derived from 
the elementary length Str^dQ)"^. The most economical way tonroceed is to compare the 
parametrization involved with that for the orthogonal ensemblecj. Two essential differences 
are easily noticed: the 2x2 blocks in matrices Yi and Y2 are interchanged and all the 
conjugated Grassmann variables have the opposite sign. The last difference, however, does 
not lead to any change in the calculation, so long as the contribution to the length Str{dQY 
from the Grassmann variables is due to terms of the kind rjdn, Tjdi] etc. Taking this into 
account, we can immediately reduce the elementary length to the following expression 



StridQY = StriidQoY + [SZ,Qo]^) , 

6Z = SR{Yo6TYo + dRR + RdSSR + 6Yo)RS , 

where all the terms apart of the last one, S R6Yq RS, are identical to those in Ref.@. 
Using Eq. (p.9|), we write 6Yq as 



(6.1) 



6Y0 = 6Y, + 6Y2 + Y1Y2SYS Y2 , 
which can be rewritten in the form 



SYn 



2 



dPiT3COs92 \ _ r ( 
dwT^w j \0T2 



(6.2) 



(6.3) 



-Ai: 



T2df32 sin 62 




d0. 



Tl 





where 



T2 




(6.4) 



Multiplying three matrices with each other, one obtains 



ro^TFn = 1x2 



■0 



cos 



+2iAi 




Tidi 
-dr]'Ti 

Tidn' \ i 

dK'Ti ; 2 



(6.5) 



where dn' = dnw exp[i(/? — /?i)/2] and dK' = wdH exp[i(/3 — /3i)/2] and analogous for dt]' 
and dr]' . One should keep in mind that the differentials dr] and dn in Eq. ( |6.5| ) are not 
the initial variables entering Eq. (|3.14| ) but new variables obtained from the initial ones 
by several replacements and shifts common for all three ensembles. The Jacobian of those 
transformations Jq is given by Eqs. ( p.lTp . 

After that we pick up the differentials of the Grassmann variables (proportional to the 
unit matrix), make a shift of the differentials analogous to the one in Ref.o and introduce 
the matrix differentials 



12 



where 



da 



da2 
da* 

dp2 
dpi 



dai da2 
—dal dal 



, dp 



dpi dp2 
-dpi dpi 



(6.6) 



— i cos ^ sinh ^drj — sin ^ cosh ^dn* 
—i cos ^ sinh ^drj* + sin ^ cosh ^dn 
—i cos Y sinh '^dn + sin y cosh ^rfr^* 



(6.7) 



-i cos Y sinh ^dn* — sin y cosh ^drj 

The Jacobian of the transformation, Eq. ( |6.7|) , from 77, k, to a and p, equals 

A 

Ju 



(cos 6*2 — cosh pY 
The supermatrix 5Z from Eq. ( |6.1|) can be represented as 
5Z = 6Y^ + iAi {2da - d9/2) + Ix 2kdp 



(6.8) 



(6.9) 



where 



1 
-1 



and 



dp sinh p 





ri 





T2 



(6.10) 



-4 



— sin 6*2 



T2 





dOo. 



n 




Calculating the length Str{dQY one notices that the anticommuting part of 5Z^ propor- 
tional to Ai, decouples from the commuting one, proportional to the unit matrix 1. The 
contribution to the length from the first part is the same as that of Ref.il leading to the 
Jacobian 



1 



1 



2^^ (sin^ + sinh^ x)^ ' 
whereas the second part of the elementary length equals 

Str[5Z\\,QQf = A{[{dpf + {d(3f sinh/x] sinh^ x + {dOf cos^ Lp 
+ {deif cosh' X + {de2f + {dpif sin' 62} . 



(6.11) 



(6.12) 



Since in our parametrization the blocks from the commuting variables in the matrices Qo and 
T are the same as in Ref.il, the Jacobian Je does not change. Combining the contributions 
from Eqs. (|6^ , |6.11| , |6.12| ) with Jg, we arrive at the elementary volume, Eq. (|3.16| - pT8|) . 
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FIGURE CAPTIONS 



Figure 1. The dependence of complex eigenvalues of energy P{e,y) on the imaginary 
part of energy y — xA/27r, for several values of non-Hermiticity a — 27rDo/i^/A. 
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